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Abstract —This paper considers three aspects of Massive 
MIMO (multiple-input multiple-output) communication net¬ 
works that have received little attention in previous works, 
but are important to understand when designing and imple¬ 
menting this promising wireless technology. First, we analyze 
how bursty data traffic behaviors affect the system. Using 
a probabilistic model for intermittent user activity, we show 
that the spectral efficiency (SE) scales gracefully with reduced 
user activity. Then, we make an analytic comparison between 
synchronous and asynchronous pilot signaling, and prove 
that the choice between these has no impact on the SE. 
Finally, we provide an analytical and numerical study of 
the SE achieved with random network deployment. 

L Introduction 

Space-division multiple access (SDMA) is a technology 
for high spectral efficiency (SE) in cellular networks, 
based on serving multiple user equipments (UEs) in the 
same time-frequency resource. While the basic concepts 
behind SDMA are more than 20 years old |[T1-||3], a 
more thorough understanding of the associated multi-user 
information theory with multi-antenna base stations (BSs) 
is a younger science ll-El. The ultimately most scalable 
form of SDMA is Massive MIMO, conceived in Q, which 
exploits the following properties Q-lil: 1) Having many 
more antennas than UEs give high spatial separability and 
near optimality of linear processing; 2) Channel state infor¬ 
mation (CSI) is only needed at the BS and can be obtained 
by uplink pilot signaling in time-division duplex (TDD) 
mode; 3) Channel hardening suppresses frequency/small 
scale fading, which makes the power allocation easier and 
eliminates the need for downlink pilots. The analytical 
study Eol has recently shown that Massive MIMO can 
provide 50-fold gains in SE over IMT-Advanced, without 
any BS cooperation, thus making it a strong candidate 
technology for next generation wireless technologies im. 

As the basic communication theoretic performance anal¬ 
ysis and asymptotic limits of Massive MIMO are becom¬ 
ing mature, it is particularly important to consider practical 
implementation aspects of the technology. The first real¬ 
time Massive MIMO testbed was described in ina and it 
shows that the Massive MIMO technology really works. 

This research has received funding from the EU FP7 under ICT-619086 
(MAMMOET), ELLIIT, and CENIIT. 


In this paper, we consider three aspects of Massive 
MIMO that are important in practical implementations; 

• Intermittent user activity: Most prior studies on 
Massive MIMO consider full-buffer transmission, 
where each UE transmits/receives data continuously. 
In contrast, practical higher layers can create bursti- 
ness (even under data streaming), thus there is not 
always data to transmit. We take a first step towards 
a more realistic traffic modeling by deriving a new SE 
expression, for a particular model of intermittent user 
activity, and analyze the performance implications. 

. Synchronous versus asynchronous pilot signaling: 
Since the seminal paper Q by Marzetta, synchronous 
pilot transmission across all cells has been the com¬ 
mon practice, in an attempt to reduce the inter-cell 
interference during pilot transmission. There are sev¬ 
eral reasons to question this basic assumption. Eirstly, 
time-synchronous transmission is hard to obtain over 
a large coverage area. Secondly, coordinated per-user 
pilot allocation is cumbersome between cells since 
the set of active UEs changes by scheduling and inter¬ 
mittent user activity, while fixed pilot reuse patterns 
with random pilot allocation within each pilot cluster 
is a more practical solution ifTOl . lfT3l . Thirdly, ifT^ 
and El showed that having unsynchronized pilot 
transmission will not cause any fundamental differ¬ 
ences. We compare synchronous and asynchronous 
pilot signaling, and obtain exciting new analytical 
results when using maximum ratio (MR) combining. 

• Performance analysis in asymmetric deployments: 
The resource allocation papers ii and El have an¬ 
alyzed the anticipated performance in symmetrically 
deployed Massive MIMO networks. In contrast, we 
analyze the SE in randomly deployed uncoordinated 
networks where the BSs are distributed according 
to a Poisson point process (PPP), which models the 
randomness of practical network deployments El- 
Closed-form SE expressions and limits are derived. 

II. System Model 

We consider the uplink of a cellular network where each 
BS is equipped with M receive antennas and serves K 
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single-antenna UEs. Since the uplink performance of a cell 
depends on the receive processing performed at the local 
BS, but is independent of the processing at other BSs, we 
can without loss of generality concentrate on BSq with 
some arbitrary cell index 0. This cell is surrounded by 
interfering cells with indices in some arbitrary set T'. 

The uplink time-frequency resources are divided into 
blocks of Tc seconds and Wc Hz, which leaves room for 
S = TcWc transmission symbols for pilot signaling and 
payload data. The channel response hijk G between 
BS/ and UE k in cell j is modeled as block fading, thus 
hijk is constant within a block and takes independent 
random realizations across blocks. More precisely, we con¬ 
sider circularly symmetric complex Gaussian realizations 
hijk ~ CA/'(0,/3/jfclM), where Pijk > 0 is the average 
channel attenuation between UE k in cell j and any of the 
antennas at BS/. This is known as Rayleigh block fading. 

Since continuous data services can give rise to intermit¬ 
tent data bursts, we assume that each UE has a probability 
A of being active in a given coherence block, where 
0 < A < 10 This is a key difference from prior works 
were only the special case A = 1 has been considered. The 
intermittent user activity of UE i in cell j is represented 
by a binary random variables aji G {0,1}, where aji = 1 
if the UE is active in the current block and aji = 0 stands 
for inactivity. These events occur with probability A and 
1 — A, respectively, and are independent between blocks. 

The received signal yo G C'^ at BSq in an arbitrary 
coherence block is then modeled as 

K K 
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where the symbol Sji € C transmitted by UE i in cell 
j has the average energy pji = E{|5jip} and no ^ 
CNiQ^a^M) denotes temporally white noise. 

III. SE UNDER Intermittent User Activity 


to channel realizations, inter-cell interference, and user 
activation variables (conditioned on aok = IJ- 

Proof: This is obtained by lower bounding the mutual 
information in the same way as in ifTOl Lemma 2], while 
taking the user activity probabilities into account. ■ 
The ergodicity of the spectral efficiency SEq/c is with 
respect to the random channel realizations and user ac¬ 
tivities. The SE depends on the combining vector got, 
which ideally should be matched to the current channel 
realizations; for example, to maximize the signal gain 
|E{hpQ^.gofe}p. Channel estimation is thus very important. 

A. Channel Estimation and Maximum Ratio Combining 
As described in Section|I] we will compare synchronous 
and asynchronous pilot transmission to gain new insights 
into the performance difference between these approaches. 
Each active UE in cell 0 sends a pilot sequence of length B 
over B instances of the system model in ([T]i, where B > K 
to allow for orthogonal sequences when all UEs are active. 
Let Sji = [sji ... G denote the combined signal 
sent by UE i in cell j during the pilot transmission in cell 
0. Eor UE k in cell 0, its signal sofe is a pilot sequence 
Sofc = ^/PaC^k, where is selected from a pilot book 

V = {vi,... ,vb} (4) 

where ||v/,|p = B, for b = 1,... ,5, and the sequences 
are orthogonal in the sense that VbiV /,2 = 0 for bif^b 2 . 

The combined received signal from the pilot transmis¬ 
sion at BSo is denoted as Yq € and given by 
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where Nq G independent CA/’(0, tT^)-entries. 

Note that an inactive UE does not send pilots nor data. 

The interfering uplink signals Sji G C^, j ^ 0, are differ¬ 
ent for synchronous and asynchronous pilot transmission. 


The considered BSq applies a linear combining vector 
go/s G C^, as yggok, to detect the signal sent by its /eth 
UE, when this UE is active (i.e., agk = !)■ The next lemma 
gives an achievable SE under intermittent user activity. 


Lemma 1. If BSq knows when its kth UE is active, a 
lower bound on its ergodic capacity is 


SEq/c — A 


^ log 2 (1 -f SINRo/s) [bit/s/Hz] (2) 


where the effective signal-to-interference-and-noise ratio 
(SINK) is given in (Oil at the top of the page and the ex¬ 
pectations E{-} and variance V{-} are taken with respect 


'with this model, the number of active UEs in a cell in a given block 
coherence has a binomial distribution: B{K,A). 


Definition 1 (Synchronous pilots). All cells send pilots 
simultaneously. In each coherence block, BSj selects K 
different pilots from V uniformly at random, such that Sji = 
Ef=i s/Ppx'ji'^b where x]i-, ■ ■ ■, Xp € {0,1} are binary 
random variables. These variables satisfy X]/)=i Xp = 1 
and x’ji = 1 occurs with probability 1/B. Since at most 
one UE per cell uses each pilot, we have xji ^ 1- 

Definition 2 (Asynchronous pilots). The interfering UEs 
transmit any random signals Sji = [sjj ... (for i = 

1,..., AT and j G tf), where Sj^ is an independent random 
variable with the average energy E{|s^jp} = Pji. 

When UE k in cell 0 is active (i.e., aok = 1), BSq can 
correlate Yq in Q with the A:th pilot sequence to obtain 




SINRofc = - 


M/3oofc 


K 


PoOk + S X) 


i=l 


ApjiPoj 

PokB 


PokB 


PoOk 


K 

E 

2=1 

i^k 


ApOi(iQQi 

POk 


K 

E E 

2=1 


Apji^Oji 

Pok 


K 


+ ^ + E E 


Ap%P^ji{M+l-A) 




( 10 ) 


K 

YoVfe =i/Pofehoofc B+Y^ YaJii^oj^s],vl + Nqv^ 

jg>i> i=i 

for fc € {!,..., R:}, (6) 

where Ngv^, ~ CA/’(0, Bct^Im)- Clearly, YqV^, is a 
sufficient statistic for estimating hoofe. Since Bayesian 
estimators require known statistics of the interfering sig¬ 
nals, which can be hard to obtain under intermittent user 
activity, we consider the least-square (LS) estimator 

hoofe = ^=%Yov:.. (7) 

yJPOkB 

Note that if UE k is inactive (i.e., a^k = 0) then hoofc = 0 
since YqV^, only contains interference and noise. 

We assume that BSq utilizes the LS channel estimate 
in 0 for MR combining where gofc = hoofc- This receive 
combining scheme attempts to maximize the desired signal 
gain |E{hQQ^gofc}p in Lemma [T] With MR we can now 
obtain closed-form expressions for the achievable SEs in 
Lemma [T] for both synchronous and asynchronous pilots. 

Theorem 1. For both synchronous and asynchronous pilot 
transmission, the lower bound on the ergodic capacity of 
UE k in cell 0 in Lemma Q] becomes 

A(l - I) log 2 (1 + SINRofe) [bit/s/Hz] (8) 

for MR combining, where the effective SINR, SINRofc, is 
given in (O at the top of the page. 

Proof: The proof is outlined in the appendix. ■ 
Interestingly, this theorem shows that the choice be¬ 
tween synchronous and asynchronous pilot transmission 
has no impact on the SE in Lemma [T] because the average 
interference during pilot transmission is suppressed by 
a factor 1/B in both cases. In the synchronous case 
this is attributed to the 1/B probability of using the 
same pilot, while the incoherent inter-cell signaling in the 
asynchronous case gives the same interference suppres¬ 
sion. Coding redundancy cannot reduce pilot contamina¬ 
tion, only increase B. The fact that also interfering data 
transmissions cause pilot contamination was previously 
demonstrated in m, ESI, thus the novel contribution of 
Theorem[T]is that the impact on the classic SE lower bound 
is exactly equal. Our numerical results hold for both cases. 

IV. Average Pereormance in Random Networks 
The performance analysis in Section|III]holds for the up¬ 
link of any Massive MIMO network deployment. To obtain 
more quantitative results, we now consider a random BS 
deployment based on stochastic geometry. More precisely, 
we assume that the BSs are distributed in according to 



Fig. 1. Example of BS positions from a homogeneous PPP $ and UEs 
uniformly distributed in the corresponding Poisson-Voronoi cells. 


a homogeneous PPP $ with intensity A BSs per km^. This 
means that in any sub-area of size C km^ there are Po(AC) 
uniformly and independently distributed BSs, where Po( ) 
is a Poisson distribution. Each UE connects to its closest 
BS, thus the coverage area of a BS is its Poisson-Voronoi 
cell; see Pig. [T] The key features of this model are the 
asymmetry and that the BSs are matched to a spatially 
heterogeneous user load, where small cells have higher 
UE density (measured per km^) than larger cells. 

The goal of this section is to characterize the SE per 
cell in the network, which can be done by studying 
the statistical properties of any UE, because each UE is 
statistically “typical” due to the translation invariance of 
homogeneous PPPs El. Using the same notation as in 
previous sections, we focus on a typical UE with index k 
that is connected to the typical BS BSq € T>. The following 
result is obtained from basic PPP properties (see El): 

Lemma 2. The distance doofc between the typical UE and 
BSq is distributed as dook ~ Rayleigh( ^ . The BSs 
of other cells form a homogenous PPP T' = $ \ {BSq} in 

\ {x e ir 2 . ii^ii ^ 

The propagation distance between BSj and UE k in cell 
j is denoted by dijk (in km). The corresponding channel 
attenuation is assumed to be distant-dependent as Pijk = 
uj~^d/j'^, where a > 2 is the pathloss exponent and oj 
determines the average channel attenuation at a reference 
distance of 1 km. Since this model causes large channel 
quality variations within the cells, uplink power control is 
necessary to avoid near-far blockage where weak signals 
from distant UEs drown in stronger signals from nearby 
UEs. We consider statistical channel inversion where UE 
i in cell j uses the average symbol energy 

Pji — (11) 

where p > 0 is a design variable and jdjji is the channel 
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attenuation from the serving BS. This power control policy 
provides the same average channel gain E = 

Mp to all UEs in the network, without having to feedback 
any instantaneous channel realizations M- The maximum 
output power of the UEs limits the choice of p in practice. 

Under these assumptions, we compute a tight lower 
bound on the average of the SE expression in Theorem [T] 

Theorem 2. A lower bound on the average uplink capacity 
[bit/symbol/user] is 

SE = a(i - log, fl + SINK ) (12) 
where SINK is given in (El at the top of the page. 

Proof: The proof is outlined in the appendix. ■ 
The closed-form average SE in Theorem|2]depends on a 
number of parameters; A, B, M, K, a, and pjcP‘. Interest¬ 
ingly, the power control policy makes the SE independent 
of the BS density A and of w, but the average energy per 
symbol depends on A as K{pji} = E{-g^} = pu!T{a/2 — 
l)/(7rA)“/^, where r(-) is the Gamma function. Most of 
these parameters are decided by the deployment scenario, 
but the number of pilot sequences, B, can be optimized. 

Corollary 1. The average uplink SE in Theorem |2] is a 
concave function of B, having its unique global maximum 
somewhere in the range K < B < S. 

Proof: The average SE has the structure A{\ — 
i3/S')log2(l -f l/(C'i + C 2 IB)), for some constants 
Cl, C 2 > 0, which is a concave function of B. ■ 

This corollary shows that there is a potential benefit 
of using pilot sequences with a length B that is longer 
than K, to suppress inter-cell interference under both 
synchronous and asynchronous pilot transmission. The 
asymptotic behavior is given by the next corollary. 

Corollary 2. The average uplink SE in Theorem\2\behaves 
as 

— ^ ■ I) + ^^AK 

when M -A 00 . This asymptotic limit is maximized by 

B = - ^ ( 15 ) 

>r(«(i + ^fea)) a-1 

where the number e is the base of natural logarithm and 
W{-) is the Lambert W function^ 

Proof: This follows from taking the limit M ^ 00 m 
(fTJI i and then finding the point where the first derivative 
of the concave SE function in (fT4l i is zero. ■ 

^The Lambert W function is defined implicitly via x = W{x)e^^A, 



Fig. 2. Empirical cumulative distribution function of the SE per cell, 
for either hexagonal and random network deployment. 


A. Numerical Examples 

We illustrate the anticipated SE in random deployments 
by a numerical study with iT = 30 UEs per cell and either 
M = 100 or M = 500 BS antennas. The asymptotic 
limit, M -A 00 , from Corollary [T] is also considered. 
The pathloss exponent is selected as a = 3.76, which 
is a common number for urban 3GPP simulations. We set 
the SNR, p/cr^, to 5 dB, which is enough to reach the 
interference-limited regime in Massive MIMO networks 
El. Each coherence block contains S = 400 symbols 
(e.g., achieved by Tc = 2 ms and Wc = 200 kHz). 

We first compare random network deployment with 
symmetric hexagonal cell deployment in terms of the SE 
per cell. Eig.|2] shows the empirical cumulative distribution 
functions (CDEs) with respect to user locations (and BS 
locations for random deployment). The curves were gen¬ 
erated by Monte-Carlo simulations based on Theorem El 
with M € {100, 500}, ^4 = 1 and B = 60, while we note 
that the average SE with random deployment is lower- 
bounded by Theorem |2] Despite the considered power- 
control, which provides all UEs with equal SNR, there are 
large SE variations since the inter-cell interference depends 
the location of the other-cell UEs. The SE variations are 
smaller with fewer antennas since the interference is more 
dominant. Hexagonal cells provide higher average SE than 
random deployment (-1-44% at M = 100 and -1-36% at 
M = 500) and gives CDEs with shorter tails, since there 
are often more neighboring cells in random deployment. 
Hence, the BS deployment has a major impact on the SE. 

Fig. [ 3 ] shows the average SE per cell, over different 
random BS and UE positions, as a function of the user 
activity A. The pilot length B is optimized numerically to 
give maximal SE. The simulation is based on the average 
of the SE in TheoremEJ Eig.[3]shows both the closed-form 
lower bound in Theorem |2] and Monte-Carlo simulations 
that lead to an upper bound due to a finite number of 
interfering cells. The small gap between the curves show 






















Fig. 3. Average SE per cell, for optimized pilot lengths, as a function 
of the intermittent user activity (i.e., the percentage of coherence blocks 
were each UE has data to transmit). There are X = 30 UEs per cell. 

the tightness of the closed-form expression in Theorem |2] 

The average SE in Fig. [3 decreases with the user 
activity, since each UE is only active in a fraction A of the 
blocks. However, the performance loss is not proportional 
to the user activity, but often much smaller. This can be 
seen by comparing the lower bound with the line “part- 
time sleep” that shows the case where only a fraction 
A of all coherence blocks are used but all K UEs are 
active in these blocks. The reason is that the network is 
interference-limited so the loss in per-cell SE from only 
having, on average, AK < K active UEs in the cell 
is partially compensated by the lower interference from 
having, on average, AK < K interfering UEs in the other 
cell. In other words, intermittent user activity seems to 
have a limited affect on the performance; for example, 
with M = 100 antennas a 30% user activity only gives a 
32% performance loss, and not 70% as one might expect. 

Adding more BS antennas improves the average SE 
in Fig. [3 The gap and relative loss as compared to the 
asymptotic limit (M — oo) is smaller for low user 
activity probabilities, which is encouraging for practical 
implementation—fewer antennas seem to be needed to 
reach a certain percentage of the asymptotic performance. 
The system achieves 30%-64% of the limit with M = 100, 
depending on A, and 64%-85% of the limit with M = 500. 

The impact of the pilot sequence length, B, is illustrated 
in Fig. |4] The average SE per cell is shown as a function 
of B, for A = 0.5 and different number of antennas. The 
optimal B is marked on each curve and, interestingly, it 
lies in the range 60-100 even if there is only K = 30 UEs 
per cell. Hence, it is preferable to have 2-3 times more 
orthogonal pilots than the total number of UEs per cell, 
even if each cell on average only needs AK = 15 of them. 
The performance variations with respect to B are quite 
small, thus any selection of B in this range 60-100 works 
well. Fig. |4] shows that the optimal B increases with M. 
The asymptotic limit (M —>■ oo) is maximized at B* = 97, 
which matches the value provided by Corollary ID 



Fig. 4. Average SE per cell, for A = 0.5, as a function of the pilot 
sequence length. There are K = 30 UEs per cell. The maximum value 
of each curve is marked. 


V. Conclusions 

We have taken a look at Massive MIMO from three 
important practical perspectives, which have been greatly 
overlooked in previous analytical works. 

First, we modeled the natural data traffic variations by a 
probabilistic intermittent user activity model, where each 
UE is on average only active in a fraction of the blocks. 
The analysis shows that reducing the user activity shrinks 
both the pre-log factor and the average interference. Since 
the latter is a positive effect, the loss in SE is smaller 
than one might expect. In fact, the simulations show that 
one needs fewer antennas to come close to the asymptotic 
limits under intermittent user activity. 

Next, we made an analytic comparison of synchronous 
and asynchronous pilot transmission and proved that it 
has no impact on the uplink performance. This confirms 
previous observations and is good news since it reduces 
the need for inter-cell pilot coordination. The result was 
proved for MR combining with random pilot allocation, 
but similar results are expected with zero-forcing pro¬ 
cessing and in the downlink. Different results might be 
obtained if well-designed pilot reuse patterns are consid¬ 
ered. Nevertheless, our results show that it is the existence 
of inter-cell interference during pilot transmission that is 
the main pilot contamination issue, while the structure of 
this interference is less important. The pilot contamination 
effect can thus be analyzed and evaluated in Massive 
MIMO implementations by sending any interfering signal. 

Finally, we analyzed the performance in randomly de¬ 
ployed networks, where the BS locations are given by a 
PER A tight lower bound on the average SE was computed 
in closed form and the limit of many BS antennas was 
characterized analytically. Simulations were provided to 
show how random deployment, intermittent user activity, 
pilot lengths, and number of antennas affect the system. 
Random deployment provides substantially lower perfor¬ 
mance than hexagonal cells, thus showing the importance 





























of coordinating the deployment of Massive MIMO BSs. we obtain an expression only containing the expectations 


Appendix A 
Outline of Proofs 

Proof of Theorem [T] The expectations in the SINR 
expression Q need to be computed for gofe = hoofc in the 
two pilot transmission cases. We begin by computing 
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by utilizing Q, the independence of the channels, and the 
fact that we conditioned on aok = 1- For the same reasons, 
the desired signal gain becomes 


|E{h«o,hoo4P = |E{h«o,hoofc}P = M^I3lok. (17) 


from which we also obtain the variance V{hQpj,hoofe} = 
E{|hoofchoofcP} - |E{h^oj.hoofc}p by noting that 

E{|h“p,hoofcp} = /3oofeE{||hoofcf} + (18) 


after some simple algebra. Next, the intra-cell interference 
terms, with i ^ k, becomes 

E{ I aoi hpQj hoofc P } = A/3ooiE{ 11 hoofc | p } (19) 


since the intra-cell pilot sequences are orthogonal. 
The inter-cell interference terms are computed as 

E{|ajihQjjhoofc|^} = A/3ojiE{||hoofe||^} 
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by keeping in mind that ho^i appears in the expression 
for hoofc. The terms E{|sJjV^p} depend on whether 
synchronous and asynchronous pilots are considered. In 
the synchronous case we have 


E{|s>:|^} = 
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since an active UE i in cell j uses the kth pilot sequence 
with probability 1/B. In the asynchronous case we notice 
that 


E{\sl^rlf}=vlE{sls*Jvl=p,.vllB^rl=p,.B, 

( 22 ) 

since Sji contains independent symbols with average en¬ 
ergy We notice from (l2Tl i and (l22l) that E{|sJjV^p} = 
PjiB in both cases, thus the SINR expression is the same. 
Finally, we obtain ( fTOl i by plugging (fT6]l- (l22l i into Q. 


Proof of Theorem |2l We compute a lower bound on 
E{SEofe}, for SEgfc given in Theorem[Tl by using Jensen’s 
inequality as 

E{SEofc} > a(i - I) log2(l + E{SINRp-i}) . (23) 

It remains to compute all expectations that appear in 
E{SINR^^^}. When applying the power-control from (fTTl i. 
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which were computed in m Proof of Prop. 1]. The 
inequality in (l25T l further establishes (fOT l as a lower bound. 
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